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ABSTRACT 

By performing a series of two-dimensional, special relativistic magnetohydro- 
dynamic (MHD) simulations, we study signatures of gravitational waves (GWs) 
in the magnetohydrodynamically-driven core-collapse supernovae. In order to 
extract the gravitational waveforms, we present a stress formula including con- 
tributions both from magnetic fields and special relativistic corrections. 

By changing the precollapse magnetic fields and initial angular momentum 
distributions parametrically, we compute twelve models. As for the microphysics, 
a realistic equation of state is employed and the neutrino cooling is taken into 
account via a multiflavor neutrino leakage scheme. With these computations, we 
find that the total GW amplitudes show a monotonic increase after bounce for 
models with a strong precollapse magnetic field (lO^^G) also with a rapid rotation 
imposed. We show that this trend stems both from the kinetic contribution of 
MHD outflows with large radial velocities and also from the magnetic contribu- 
tion dominated by the toroidal magnetic fields that predominantly trigger MHD 
explosions. For models with weaker initial magnetic fields, the total GW ampli- 
tudes after bounce stay almost zero, because the contribution from the magnetic 
fields cancels with the one from the hydrodynamic counterpart. These features 
can be clearly understood with a careful analysis on the explosion dynamics. 
We point out that the GW signals with the increasing trend, possibly visible to 
the next-generation detectors for a Galactic supernova, would be associated with 
MHD explosions with the explosion energies exceeding 10^^ erg. 

Subject headings: supernovae: collapse — gravitational waves — neutrinos — 
hydrodynamics 
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Introduction 



Successful detection of neutrinos from SN1987A paved the way for Neutrino Astronomy 
( iHirata et al.lll987l : iBionta et al.lll987l ) , alternative to conventional astronomy by electromag- 
netic waves. Core-collapse supernovae are now expected to be opening yet another astron- 
omy, Gravitational- Wave Astronomy. Currently lon g-baseline laser interferometers LIGO 
jAbbott et aljboosh . VIRGO\ GEOeOO^, TAMA300 (lAndo fc the TAMA Collaboration|[2005l ) 
and AIGO^ with their international network of the observatories, are beginning t o tak e data 
at sensitivities where astrophysical events are predicted (see, e.g., iHough et al.l (120051 ) for a 
recent review). For these detectors, core-collapse supernovae have been proposed as one of 
themost plausible sources of gravitational waves (GWs) 
(120091 ) for recent reviews). 



(^see, e.g. 



Kotake et al.l (120061 ): lOtt 



Although the explosion mechanism of core-collapse supernovae has not been completely 
clarified yet, current multi-dimensional simulations based on refined numerical models show 
several promising scenarios. Among the candidates is the neutrino hea ting mechanism 
aided b y convection and sta n ding accretion shock insta b ility ('SASI) (e.g., iMarek fc Janka 



nism (e.g.. 



mechanis m (iBurrows et al. 



mm -. iBruenn et al.l (120101): IScheck et al.l (12004 |2008|); ISuwa et al.l (l2009l l). the acoustic 



20061 . l2007br). and the magneto h ydrodynamic (MHD) mecha 



Ardelian et al.l (120001): iKotake et al.l ( l2004aL l2005l l: lObergaulinger et al.l ( ]2006br ): 



Burrows et all J2007ah : iTakiwaki et al.l ( l2009l ) and references therein). For the former two 
to be the case, the explosion geometry is expected to be unipolar and bipolar, and for the 
MHD mechanism to be bipolar. 

Since the GW signatures imprint a live information of the asphericity at the moment of 
explosion, they are expected to provide us an important hint to solve the supernova mecha- 



in the context of rotational core-collapse (e.g.. 


Monchmever et al. (1991): 2 


jwereer & 


Mueller 


(1997): 


Kotal 


le et al. 


(2003b): Shibata k Sekieuchi (2004): Ott et al. 


(2004 


. 2007a, b): 


Dimmelmeier et al. 


(2002a, 


2007, 


2008): 


Scheideeeer et al. (20081 


). For the bounce signals having a strong and 



characteristic signature, the iron core must rotate enough rapidly. The waveforms are cate- 
gorized into the following three types, namely types I, II, and III. Type II and III waveforms 
are shown less likely to appear than type I, because a combination of general relativity (GR) 
and electron capture near core bounce suppresses multiple bounce in the type II waveforms 



http://www.ego-gw.it/ 
^http://geo600. aei.mpg.de/ 
•^http:/ /www. gravity, uwa.edu.au/ 
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( iDimmelmeier et al.l 120071 : lOtt et al.l l2007al Jbl). In general, a realistic nuclear equation of 
state (EOS) is stiff enough to forbid the type III waveforms. After bouiice, asymmetries du e 



to con vection (IBurrows &: Hayes 



SASI flMarek et al 



2009 



199 



i 



Kotake et al 



200 



2009b 



a 



Miller &: Janka 

i 



1997 



Frvei]l2004j: Miiller et al.ll2004[) 



Murphy et al. 20091 ) . and g-mode os- 



cillations of protoneutron stars ( Ott et al. 20061 ). are expected to account for sizable GW 
signals. 

In general, detection of these GW signals in the postbounce phase (except for the g-mode 
oscillation) is far more difficult than the bounce signals, because they do not possess a clear 
signature like bounce signals, but change stochastically with time as a result of chaotically 



grow i ng convection as wel l as SASI in the non-l inear hydrodynamics (iKotake et al.l ( l2009b 
201lh : iMarek et all J2OO9I ): iMurphv eTaP toO^ ). 



Rapid rotation, necessary for the strong bo u nce s ignals, is likely to obtain ~ 1% of 
massive star population (e.g., IWoosley fc BloomI ( l2006l )). However this can be really the 
case for progenitors of rapidly rot ating metal-poor stars , which experience th e so-called 
chemically homogeneous evolution (IWoosley fc Hegerll2006l : lYoon fc Lange ] l2005h . The high 
angular momentum of the core as well as a strong precollapse magnetic field is preconditioned 
for the MHD mechanism, because the MHD mechanism relies on the extraction of rotational 
free energy of the collapsing core via magnetic fields. The energetic MHD explosions are 
receiving grea t atten t ion re cently as a possible relevance to magnetars and collapsars (e.g., 
Harikae et al.l ( 120091 . |2010|) for collective references), whic h are pres u mably linked to the 



formation of long-duration gamma-ray bursts (GRBs) (e.g., iMeszarosI ( 120061 )) 



Among the previous studies mentioned above, only a small po rtion of papers ha s been 



Kotake et al.l 


200^ 


lb; 


Obereaulineer et al. 


2006b 


Scheidegger et al. 


I2OI0I). This may be because 



MHD mechanism 


(Yamada & Sawai 




2004; 


Cerda-Duran et al. 


2007 


Shibata et al. 


2006 



their infiuence over the GW si gnals can be visible only for cores with precoll apse mag- 
netic fields over Bq > 10^^ G ( lObergaulinger et al.1 l2006bl : Kotake et al.ll2004b[ ). Consid- 



eri ng that the typical mag netic-field strength of GRB progenitors is at most 



10 



11-12 



G (IWoosley &: Hegerl |2006| ). this is already an extreme situation. Interestingly in a more 

Q-v-fvcin-i £il T r o 7-?_ « , 1 O-^^ n o/Ui^n 1 o vIt r rfrrstnTTr^rf f o n v/ii m ■J-Vi/Zi Tirm r/zifr^vwi o Tiro o r\\~\ 



served ( 


— "<i 
Obereaulineer et al. 


„ , 
2006a; 


Obergaulinger et al. 


( 


2006a 


called 



scale oscillations of GWs near bounce are replaced by higher frequency irregular oscillations. 
Some of these MHD simulations follow adiabatic core- collapse, in which a polytropic EOS 
is employed to mimic supernova microphysics. At this level of approximation, the bounce 
shock generally does not stall and a prompt explosion occurs within a few ten milliseconds 
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after bounce. Therefore a main focus in these previous studies has been rather hmited to 
the early postbounce phase (< several 10 ms). However, for models with weaker precol- 
lapse magnetic fields akin to the current GRB progenitors, the prompt shocks stall firstly 
in th e core like a conyentiori al supernova model with more sophisticated neutrino treatment 
(e.g., iBurrows et al.l ( l2007al )). In such a case, the onset of MHD exp losions, depending on 
the initial rotation r ates, can be delayed till ~ 100 ms after bounce ( iBurrows et al.ll2007al : 



Takiwaki et al.ll2009l ). There remains a room to study GW signatures in such a case, which 
we hope to study in this work. 

In this study, we choose to take precollapse magnetic fields less than 10^^ G based on a re- 
cent GRB-oriented progenitor models. By this choice, it generally takes much longer time af- 
ter bounce than the adiabatic MHD models to amplify magnetic fields enough strong to over- 
whelm the ram pressure of the accreting matter, leading to the magnetohydrodynamically- 
driven (MHD, in short) explosions. Even if the speed of jets in MHD explosions is only mildly 
relativistic, Newtonian simulations are not numerically stable because the Alfven velocity 
(oc B/ y/p) could exceed the speed of light unphysically especially when the strongly magne- 
tized jets (i.e., large B) propagate to a stellar envelope with decreasing density (p). To follow 
a long-term postbounc e evolution nurnerical ly stably, we perform special relativistic MHD 
(SRMHD) simulations (ITakiwaki et al.ll2009l ). in which a realistic EOS is employed and the 
neutrino cooling is taken into account via a multiflavor neutrin o leakage scheme. N ote in 
our previous study of GWs in magneto-rotational core-collapse (IKotake et al.ll2004bl ) that 
we were unable to study properties of the GWs long in the postbounce phase because the 
employed Newtonian simulations quite often crashed especially in the ca se of strong MHD ex 



plosion s. To include GR effects in this study, we follow a prescription in lObergaulinger et al. 



(j2006al ) which is reported to capture basic features of full GR simulations quite well. By 
changing precollapse magnetic fields as well as initial angular momentum distributions para- 
metrically, we compute twelve models. By doing so, we hope to study the properties of GWs 
in MHD explosions systematically and also address their detectability. 

The paper opens up with descriptions of the initial models and numerical methods 
employed in this work (section [2]) . Formalism for calculating the gravitational waveforms in 
SRMHD is summarized in section |3l The main results are given in Section HI We summarize 
our results and discuss their implications in Section El 
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2. Models and Numerical Methods 
2.1. Initial Models 

We make precollapse models by taking the profiles of density, inte rnal energy, and 



electro n fraction distribution from a rotating presupernova model of E25 in iHeger fc Langer 



(|2000| ). This model has mass of 25Mq at the zero-age main sequence, however loses the 
hydrogen envelope and becomes a Wolf-Rayet (WR) star of 5.45 M© before core-collapse. 
Our computational domain involves the whole iron-core of l.GQM©. Note that this model is 
suggested as a candidate progenitor of long-duration GRBs because type Ib/c core-collapse 
supernovae originated from WR star s have a observational association with long-duration 



GRBs (e.g., IWooslev fc BloomI (120061 ) I 



Since little is known about the spatial distributions of rotation and magnetic fields in 
evolved massive stars, we add the following profiles in a parametric manner to the non- 
rotating core mentioned above. For the rotation profile, we assume a cylindrical rotation 
of 

where Q is angular velocity and X and Z denotes distance from the rotational axis and 
the equatorial plane, respectively. The parameter Xq represents the degree of differential 
rotation, which we choose to change in the following three ways, 100km (strongly differ- 
ential rotation), 500km, (moderately differential rotation), and 2000km (uniform rotation), 
respectively. The parameter Zq is fixed to 1000km. 

Regarding the precollapse magnetic field, we assume that the magnetic field is nearly 
uniform and parallel to the rotational axis in the core and dipolar outside. This can be 
modeled by the following effective vector potential, 

Ar = Ae = 0, (2) 

A, = ^^jsin9, (3) 
where Ar,e,</. is vector potential in the r, 9, cj ) direction, resp e ctivel y, r is radius, tq is radius 



of the core, and Bq is a model constant (see iTakiwaki et al.l (120041 ) for detail). In this study. 



ro is set to 2000 km which is approximately the size of the precollapse iron core. 

By changing initial angular momentum, degree of differential rotation, and the strength 
of magnetic fields, we compute twelve models. The model parameters are shown in Table [H 



- 6 - 



The models are named after this combination, with the first letters, B12, Bll representing 
strength of the initial magnetic field, the second letters, X1,X5,X20 indicating the degree 
of differential rotation {Xq = 100, 500, 2000 km, respectively), and the third letter, /3 = 0.1, 1 
showing the rotation parameter /3. Here /3 represents ratio of the rotational energy to the 
absolute value of the gravitationa l energy prior to core-collapse. The original progenitor 
of model E25 in iHeger &: Langed (120001 ) has a uniform rotation profile in the iron core 
and the initial (3 parameter is ~ 0.15 %. So the initial angular momentum in our models of 
B11X20/30.1 and B12X20/30.1 are sirn ilar to the original one. In a GRB-oriented progenitor of 
model 350B (IWoosley &: Hegerll2006l ). the precollapse magnetic fields reach to ~ 10^^ — 10^^ 
G and /3 ~ 0.2%, which is not so different from the chosen parameters here. Although little 
is known about the degree of precollapse differential rotation, an extremely strong one, for 
example, our model series of XI, should be unrealistic. These models, albeit rather academic, 
are examined in order to see clearly the effects of differential rotation. 





100km 


0.1% 
500km 


2000km 


/3(%) 
Xo(km) 


100km 


1% 

500km 


2000km 


Bo : lO^^G 
lO^^G 


BllXl/30.1 
B12X1/30.1 


B11X5/30.1 
B12X5/30.1 


B11X20/30.1 
B12X20/30.1 




BllXl/31 
B12X1/31 


B11X5/31 
B12X5/31 


B11X20/31 
B12X20/31 



Table 1: Summary of initial models. Model names are labeled by the precollapse magnetic 
fields and rotation. /3 represents ratio of initial rotational energy to the absolute value of the 
initial gravitational energy. From left to right in the table, Qq in unit of rad/s (equation ([1])) 
is 24, 2.8, 0.95, 76, 8.9, and 3.0, respectively. Note that Xq and Bq is defined in equation (1) 
and (I3j), respectively. 



2.2. Special Relativistic Magnetohydrodynamics 



Numerical results in this work are calculated by the SRMHD code developed in lTakiwaki et al 



(120091 ). In the following, we first mention the importance of SR and then briefly summarize 



the numerical schemes. 

The Alfven velocity of MHD jets propagating into the outer layer of the iron core can 
be estimated as va ~ 10"'^''cm/s (i?/10"'^^G)(p/10^g/cm'^) ^^"^ , where p and B are the typical 
density and the magnetic field near along the rotational axis. It can be readily inferred that 
the Alfven velocity could exceed the speed of light unphysically in Newtonian simulations. 
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SR corrections are also helpful to capture correctly the dynamics of infalling material in the 
vicinity of the protoneutron star, because their free-fall velocities and rotational velocities 
become close to the speed of light. Such conditions are quite ubiquitous in MHD explosions. 
Even if the propagation speeds of the jets are only mildly relativistic, we have learned that 
(at least) SR treatments are quite important for keeping the stable numerical cal culations 
in good accuracy over a long-term postbounce evolution (e.g., iHarikae et al.l (120091 )). 



The MHD part of our code is based on the formalism of iDe Villiers et al.l ( 120031 ). The 
state of the relativistic fluid element at each point in the space time is described by its 
density, p; specific energy, e; velocity, f*; and pressure, p. And the magnetic field in the 
laboratory frame is described by the four-vector V^6^ = *F^'^Uy^ where *F^'^ is the dual of 
the electro-magnetic field strength tensor and is the four-velocity. The basic equations of 
the SRMHD code are written as, 



9D 1 ^ i 

a pi 1 

d{S,-h%) , 1 . 

dt + ^^.v^(^^-^-^^^) 



d{WV - Wb^v'^ 
dt 



-p- 



dW 



dt 



(4) 
(5) 



--{ph{Wvj.f 



kk 



-di{p + 



+ dj{Wv^U -WvV) = 



(6) 
(7) 



An 



{{by + (bkr) 



where W 



D = pW, E = eW and Si = phW^Vi are the Lorentz boost factor. 



auxiliary variables correspond to density, energy, and momentum, respectively. All of them 
are defined in the laboratory frame. Equations (Hf5|6|) represents the mass, energy, and 



momentum conservation, respectively. In equation (jH]), note that the relativistic enthalpy, 
h = {l + e/p+p/p+ \bf / p) includes the magnetic energy. Equation ([7]) is induction equation 
in SR. In solving the equation, the meth od of characteristics i s implemented to propagate 
accurately all modes of MHD waves (see iTakiwaki et al.l ( 120091 ) for more detail). Equation 
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(IH]) is Poisson equation for the gravitational potential of $, which is solved by the modified 
incomplete Cholesky conjugate gradient method. For an approximate treatmei it of general 
relatiyistic (GR) gravity, $tot in equation includes a GR correction to $ as in lBuras et al. 

We empl oy a realistic equation of state based on the relativistic mean field theory 
jshen et alilim . 



We approximate the neutrino co oling by a neutrino multiflavor leakage scheme (lEpstein fc Pethick 



198ll : iRosswog fc Liebendorferl 120031 ). in which three neutrino flavors: electron neutrino (z/e), 
electron antineutrino(z/e), and heavy-lepton neutrinos (i/^, z/^, z/^, Ur, collectively referred 



to as ux), are taken into account. The implemented neutrino reactions are electron cap- 



Takahashi et al. 


1978; 


Itoh et al. 


1989. 


1990) 



the lepton fractions (namely — Y"p+ , Yj^, , Yp , and , are solved in an operator-splitting 



manner (see equation (7) in iTakiwaki et al.l (120091 ) for more detail). in equation ([5]) 
represents the neutrino cooling r ate summed over a l l the r eactions, which can be also est i- 
mated by the l eakage scheme (see Epstein fc Pethickl ( 1l98ll ) ; iRosswog fc Liebendorferl ( l2003l ) ; 
Kotake et al.l (l2003al ) for more detail). 



In our two dimensional simulations, the spherical coordinates are employed with 300(r) 
X 60{6) grid points to cover the computational domain assuming axial and equatorial sym- 
metry. The radial grid is nonuniform, extending from to 4000km with finer grid near the 
center. The finest grid for the radial direction is set to 1km. The polar grid uniformly 
covers from 6' = to 7r/2. The finest grid for the polar direction is 25 m. The numerical 
tests and conve r gence with this choice of the numerical grid points are given in section 5 of 
Takiwaki et"aD toO^ . 



To measure the strength of explosion, we define the explosion energy as follows. 



exp 



dV ei 



ocal 



D 



dV (Ckin + Cint + Cmag + 6, 



grav; 5 



(9) 



here eiocai is the sum of Ckin, Cint, Cmag and Cgrav with being kinetic, internal, magnetic, and 
gravitational energy, respectively defined as 



ekin = pW{W-l), 
eint = eW^+p{W^-l) 
1 



|6|M1 



^02 



2Vr2/ 2W^'' 



--grav 



(10) 
(11) 

(12) 

(13) 



and D in equation (9) represents the domain where the local energy (eiocai) is positive. 
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indicating that the matter is gravitationally unbound. The explosion energy is evaluated 
when the MHD jets pass through 1000 km along the polar direction. 



Formulae for Gravitational Waves in SRMHD 



To extract the gravitational waveforms in MHD explosions, we present the stress formula 
in SRMHD for later convenience. As shown below, this can be done s traight forwardly by 
extending the Newtonian MHD formulation presented in lKotake et al.l (l2004bl ). 



From the Einstein equation, one obtains the following formula as a primary expression 
for the leading part of the gravitational q uadrupole field emit t ed by t he motion of a. fluid i n 
the post-Newtonian a pproximation (e.g., iMonchmeyer et al.l ( 1l99ll ): iFinn fc Evand ( 1l990l ): 
Blanchet et aP Jl99ol )). 

AG 



c^R 



■kh 



(14) 



where G and c are the gravitational constant and the velocity of light, respectively, Tki 
is the energy momentum tensor of the source, R = {dijX'X^y/'^ = \X\ is the distance 
between the observer and the source. Pijki, with N = X / R denotes the transverse-traceless 
(TT) projection operator onto the plane orthogonal to the outgoing wave direction N (e.g., 
Monchmeyer et al.l (1l99ll )) . 



Tij consists of the three parts, namely of perfect fluid, electromagnetic fleld, and gravi- 



tational potential as follows. 



The flrst term which we refer as the hydrodynamic part is explicitly written as. 



^»i(hyd) 

where is effective density deflned as, 

P* = 



P + 



e + p+ 



(15) 



(16) 



(17) 



The second term in equation (15) represents the contribution from magnetic flelds as. 



T- ■ 

*J (mag) 



-bibj. 



And the last term in equation (15) is the contribution from the gravitational potential, 

1 . 1 



T- 

■J- 7 



^3 (grav) 



(19) 
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where $ corresponds to the self-gravity in equation ([8]). 



In our axisymmetric case, there remains only one non- vanishing quadrupole term in the 
metric perturbation, namely £ = 2,m = in terms of the pure-spin tensor harmonics, as 



TTi V 4.\ «=^_>!"=^ —A^'^ [t— — \ T-^^'^° 



R 



{OA), 



where T. 



E2,20/ 
ij \ 



IS 



pE2,20j 
-ij \ 



1 /15 




sin^ e. 



TX 



e.g., iThornd (119801 )). The projection operator in equation ( IT^ acts on Tij as. 



7~) rpkl OT" 



ZZ rrixx 



(20) 



(21) 



(22) 



Transforming equation f|T^ to the spherical coordinates, and expressing 6j and f , in terms 
of unit vectors in the r, 6, direction, we obtain for the expression. 



where 



aE2 

^20 (hyd) 



E2 

20 (hyd) 



aE2 _ aE2 , 4E2 , /(E2 

^20 ~ ^20 (hyd) ^ ^20 (mag) ^20 (grav) ' 



G 327r3/2 ^1 



JO 



20 (hyd)' 



p*W^'(t;.'(3^2 - 1) + t;e'(2 - 3^') - v^^ - 6t;,t;e /i^l-Zi^); 



(23) 

(24) 

(25) 



and 



A 



E2 

20 (grav) 



/E2 



20 (grav) 



A 



c 



~ VTE Jo 

ph{W^ + {Vklcf) + I + ^ ) - ^ + {hf) 



20 (grav) ' 
2' 

2 



X 



-(9^$(3/i2 - 1) + ^de^ 



E2 

20 (mag) 



G 327r3/2 f\ r 2, , 
djj, f r drj. 



fE2 
J 2 



rA 

2io ..2 



%/15 JO JO 



E2 

20 (mag)' 



'20 (mag) - [6,^(3^^ - 1) + he\2 - Sfi') - - 6KbefxVl - /i^] 
where fi = cos 9. For later convenience, we write the total GW amplitude as, 

_ , TT , . TT , . TT 

~ "'(hyd) "T "-(mag) "r "-(grav)' 



(26) 



(27) 

(28) 
(29) 

(30) 



- 11 - 



where the quantities of the right hand are defined by combining equations fl20|) and ( !23|) with 
equations fl25|) . fl27j) . and f l29|) . By dropping (v/c) terms, the ab o ve for mulae reduce to the 
conventional Newtonian stress formula (e.g., iMonchmeyer et al.l (Il99ll )). In the following 
computations, the observer is assumed to be located in the equatorial plane {6 = 7r/2 in 
equation (12T!) ). and that the source is located near at our galactic center {R = 10 kpc). 
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Results 



The gravitational waveforms obtained in this work can be categorized into two, which 
we call increasing type or cancellation type just for convenience. Note that the latter type 
does not mean a new waveform as will be explained later in t his section. Regarding the 
former one, such a waveform was presented i n previous literature (lObergaulinger et al.ll2006al : 



Shibata et al.ll2006l : IScheidegger et al.ll2010l ). however their properties have not been clearly 



understood yet. In section WA] we first overview their characters, which are peculiar in the 
case of MHD explosions. In section 14. 2[ we move on to analyze their properties by carefully 
comparing each contribution in equation (15U]) to the total GW amplitudes. Then in section 
14. 3[ we perform the spectra analysis and discuss their detectability. 



4.1. Properties of Waveforms in the MHD Exploding Models 

Figure [1] shows examples of the two categories, which we call as the increasing (left 
panels) or cancellation type (right panels), respectively. In the increasing type, the total 
wave amplitudes (red line) have a monotonically increase trend after bounce {t — tb = 
in the figures). While in the cancellation type (right panels), the total amplitudes after 
bounce stay almost zero. This is because the contribution from the magnetic fields (blue line, 
equation fl29|) ) cancels with the one from the sum of the hydrodynamic and gravitational parts 
(green line, equations fl25|27p ). Regardless of the difference in the two types, it is common 
that the magnetic contribution (blue line) increases almost monotonically with time. Not 
surprisingly, the bounce GW signals ((t — tf, < 20 ms) are categorized into the so-called type 
1 or 11 waveforms. Note here that the MHD simulations are terminated at around 100 ms 
after bounce for all the computed models. This is simply because the GW amplitudes in a 
more later phase decrease because the MHD shock comes out of the computational domain 
and the enclosed mass in the domain becomes smaller. 

Table [2] depicts a classification of the computed models, in which "C" and "I" indicates 
the cancellation and increasing type, respectively. 1* in the table indicates the mixture of the 
two types, which we call as intermediate type. The table shows that the bifurcation of the 
two types is predominantly determined by the precollapse magnetic fields, so that the models 
with stronger magnetic fields {Bq 10^^ G) are basically classified to the increasing type. For 
models colored by orange, which are slow rotator with uniform rotation in our models, the 
field amplification works less efficie ntly than for mode l s with stronger differential rotation 



(such as Xq = 100, 500 km) (see also lCerda-Duran et al.l (120081 )). This suppresses the increase 



in the wave amplitudes due to magnetic fields, which gives rise to the intermediate state 
between the two types. 
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Fig. 1. — Gravitational waveforms with the increasing (left) or the cancellation trend (right) 
(see text for more detail). At the right bottom in each panel, the model names are given 
such as B12X1/30.1 (top left, for example). The total wave amplitudes are shown by the red 
line, while the contribution from the magnetic fields and from the sum of the hydrodynamic 
and gravitational parts are shown by blue and green lines, respectively (e.g., equation (|29ll 
and equations ( 125|271) ). 
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Table |3] is a summary showing intervals measured from the stall of the bounce shock 
(top) and from core bounce (b ottom) till the MHD -driven revival of the stalled bounce 



shock. As already mentioned in iTakiwaki et al.l ( 120091 ). generation of the postbounce MHD 



jets proceeds in the following two ways. One is launched relatively promptly after the stall 
of the bounce shock, typically earlier than ~ 30 ms (see models colored by red in Table 3) 
and another is launched rather later after bounce (> 30 ms) (models colored by green in 
Table 3). In this sense, our computed models could be roughly categorized into two, namely 
promptly MHD explosion (colored by red in Table 3) or delayed MHD explosion (colored by 
green in Table 3), respectively. This simply reflects that it takes longer time for the weakly 
magnetized models to amplify the magnetic pressure behind the stalled shock enough strong 
to overwhelm the ram pressure of accreting matter. By comparing Table [2] to [3l the two 
characters in the waveforms have a rough correlation with the difference of the explosion 
dynamics. 

TableHlshows a summary regarding the explosion energy (e.g., equation Q)). Comparing 
Table [2] to m the explosion energy for the increasing type (models colored by yellow in Table 
2 and 4) is higher compared to the cancellation type (models colored by light blue). Given 
the rotation rate of /3 = 0.1%, the explosion energy is smallest for the intermediate type 
(models colored by orange) due to the insufficient field amplification as mentioned above. 
It can be also shown that the explosion energies for the models with the increasing trend 
generally exceed 10^^ erg (models colored by yellow). 

Having summarized the waveform classification together with the explosion dynamics, 
we move on to look more in detail what makes the difference between the two types in the 
next section. 
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Table 2: Same as Table [H but for the classification of the computed models. "C" and "I" 
indicates the cancellation and increasing type, respectively, while I* indicates the mixture of 
the two types, which we refer to as intermediate type. 
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Table 3: Same as Table [T] but for the interval measured from the stall of the bounce shock 
to the MHD-driven revival of the stalled shock (top panel) and the one measured from core 
bounce (bottom panel). The computed models are classified whether the launch of the MHD 
jets occurs relatively promptly after bounce (models colored by red, with the intervals being 
shorter than ~ 30 ms typically) or rather later (models colored by green), which we refer 
to as promptly or delayed MHD explosion for convenience in this work (see text for more 
detail). Note that these timescales are estimated just by looking at the velocity evolutions 
along the polar axis. 
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Table 4: Same as Tabled] but for the explosion energy (defined in equation (9)). Comparing 
with Table [2] to HI the explosion energies for the models with the increasing trend generally 
exceed 10^^ erg (models colored by yellow), i 
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4.2. Analysis of Waveforms 

4.2.1. Increasing type 

By taking model B12X1/90.1 as a reference, we first focus on the increase- type waveform. 
The left panel of Figure |2] depicts distributions of entropy (left-half) and plasma /3 (right- 
half, the ratio of matter to magnetic pressure) at 100 ms after bounce. It can be seen that 
the outgoing jets indicated by the velocity fields (arrows in the left-half) are driven by the 
magnetic pressure behind the shock (see bluish region (i.e., low plasma /3) in the right-half 
panel) . 

The right panel of Figure [2] shows contributions to the total GW amplitudes (equa- 
tion (27)), in which the left-hand-side panels are for the sum of the hydrodynamic and 
gravitational part, namely log ^± /^^(^yd) + /^^(grav) ) (^^^^ top(+)/bottom(—) (equations 

( I25f27|) ). and the right-hand-side panels are for the magnetic part, namely log ^±/^^(mag)) 
(right top(+)/bottom(— )) (e.g., equation ( 129|) ). By comparing the top two panels in Fig- 
ure 2, it can be seen that the positive contribution is overlapped with the regions where 
the MHD outfiows exist. The major positive contribution is from the kinetic term of the 
MHD outfiows with large radial velocities (e.g., +p^.W'^Vr'^ in equation (l25ll ). The magnetic 
part also contributes to the positive trend (see top right-half in the right panel (labeled by 
mag(-|-))). This comes from the toroidal magnetic fields (e.g., in equation f l29|) ). which 
dominantly c ontribute to drive M HD explosions. The magnetic contribution was already 



mentioned in lKotake et al.l f 2004bl ). This study furthermore adds that the kinetic energy of 



MHD outfiows more importantly contributes to the positive trend. 

Figure [3] shows a normalized cumulative contribution of each term in , which is 
estimated by the volume integral of A^q within a given sphere enclosed by certain radius. It 
can be seen that the contribution of the hydrodynamic and gravity parts (indicated by "hyd 
& grav") is prominent for radius outside ~ 1000 km, which stems from exploding regions 
with large kinetic energy as mentioned above. Note that the secular drift observed in the 
increase-type waveform may come fro m ambiguities in est i mating the gravitational potential 



in the stress formula as pointed out by lDimmelmeier et al.l ( l2002bl ). To exclude such potential 



issues , we plot the w a yeform calculated by the first-moment of momentum-density formalism 



(e.g.. 



Finn fc EvansI (Il990[ ) ) . As is shown in Figure HI the increasing trend is also seen in 
the waveform estimated by the first-moment formalism (green line). As a side- remark, a 
more smoother curve is obtained for the stress formula (red line), about which Monchmeyer 
et al. (1991) pioneeringly mentioned that the numerical evaluation of the time-derivative 
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sometimes makes the waveform noisyll] 



4-2.2. Cancellation type 

Now we proceed to focus on the cancellation-type waveform by taking model BllXl/91 
as an example. Top panel in Figure [S] shows density distributions (left-half) with the plasma 
/3 (right-half). Note that the jet head of MHD outflows is located at ~ 400 km along the 
pole which sticks out of the plots in Figure [5] (compare the difference in scales of Figure |2]). 
In fact, the bluish regions (low /3) around the pole have a lower density because material 
there has been already blown up due to the passage of MHD-driven shocks. The middle 
panel of Figure 5 shows the term-by-term contribution to A^q as in the right panel of Figure 
[21 From the right-half panel, it is shown that the magnetic contribution dominantly makes 
a positive contribution (labeled by mag(-|-)), which is also the case of the increasing type as 
mentioned in the previous section. Note that the negative contribution from magnetic fields 
(bluish region in the bottom-right-half of the middle panel) comes from the regions, where 
are off-axis from the propagation of the MHD jets. In these regions, the poloidal components 
of magnetic fields are stronger than the toroidal ones, which makes the negative contribution 
mainly through -[^^^(S/i^ - 1) + 6/(2 - S/i^)] in equation (|29l). 

Looking at the sum of the hydrodynamic and gravitational part (left-half in the middle 
panel), a large negative contribution comes from regions near in the rotational axis (colored 
by red, bottom-left-half). The bottom panel of Figure [5] further shows a contribution from 
the hydrodynamic (left-half) and gravitational part (right-half), separately. Regarding the 
hydrodynamic part, the negative contribution is highest in the vicinity of the equatorial 
plane which closely coincides with the oblately deformed protoneutron star (colored by red, 
the bottom- left-half in the bottom panel). This is because the negative contribution comes 
from the centrifugal forces (e.g., the term related to the rotational energy, —pv^"^ in equation 
f l25|) ) . For the gravity part (right-half in the bottom panel) , a big negative contribution comes 
from regions in the vicinity of the rotational axis (bottom-right-half). This comes from the 
term of — r9r$(3/i^ — 1) in equation (125|) . which is the radial gradient of the gravitational 
potential. Remembering that > is generally satisfied in self-gravitating objects, the 
sign of this term is determined by /i = cos 6 which is a directional cosine measured from the 
rotational axis. As a result, the gravity part makes a negatively contribution in the vicinity 
of the rotational axis (/i ~ 1). These two factors coming both from the hydrodynamic and 
gravity part make the negative contribution to the GW amplitudes. In fact, the negative 



This may be the reason why the stress formula has been often employed in supernova researches so far. 
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Fig. 2. — Left panel shows distributions of entropy [fc^/baryon] (left) and logarithm of 
plasma /3 (right) for model B12X1/30.1 at 100ms after bounce. The white arrows in the left- 
hand side show the velocity fields, which are normalized by the scale in the middle left edge 
(0.5c). Right panel shows the sum of the hydrodynamic and gravitational parts (indicated 
by "hyd and grav" in the left-hand side) and the magnetic part (indicated by "mag" in the 
right-hand side), respectively. The top and bottom panels represent the positive and negative 
contribution (indicated by (+) or (-)) to Afg, respectively (see text for more detail). The 
side length of each plot is 4000(km)x8000(km). 
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Fig. 3. — Normalized cumulative contribution of each term in Afg as a function of radius for 
model B12X1/30.1. This is estimated by the volume integral of Afg within sphere of a given 
radius. 
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Fig. 4. — Gravitational waveform extracted either by the stress formula (indicated by QPM, 
red line) or by the first-moment of momentum-density formalism (FDM, green line) for model 
B12X1/30.1), respectively. 
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contribution is highest near at the surface of the protoneutron star (~ 10 km in Figure |6]). 
Outside it, the magnetic part becomes almost comparable to the sum of the hydrodynamic 
and gravity part, which makes the cancellation type shown in the right panels of Figured] 
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Fig. 5. — The top panel shows distributions of logarithm of density [g/cm^] (left-half) and 
logarithm of /3 (right-half) for model BllXl/31 at 100 ms after bounce. Like Figure 121 the 
middle panel shows the sum of the hydrodynamic & gravitational part (left-half), while 
the magnetic contribution is shown in the right-half panel. The bottom panel is for the 
hydrodynamic (left-half) and gravity part (right-half), respectively. The side length of each 
plot is 400(km)x400(km). Compared to Figure 2, note that a more central region is focused 
here because the contribution to the GWs there are more important for the cancellation-type 
waveform. 
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Fig. 6. — Same as Figure [3] but for model BllXl/31. 
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4.3. Spectrum Analysis 

Now we move on to perform a spectral analysis (Figure [7]). Both for the two types (left 
panel (increasing-type), right panel (cancellation-type)), the peak amplitudes in the spectra 
are around 1 kHz, which comes from the GWs near at bounce. It can be seen that the 
spectra for lower frequency domains (below ~ 100 Hz) are much larger for the increasing 
type (left panels) compared to the cancellation type (right panels). This reflects a slower 
temporal variation of the secular drift inherent to the increase- type waveforms (e.g.. Figure 

As a measure to characterize the dominance in the lower frequency domains, we define 
hiovj, which represents average amplitudes below 100 Hz (see Table [5]). Although the peak 
amplitudes, /ipeak, in the spectra have no clear correlation with the two types, we point 
out that the final GW amplitudes (the first column) and the /iiow (the third column) are 
much larger for the increasing type (colored by yellow) compared to the cancellation type 
(colored by light blue). In Figure 7, the peak amplitudes near 1 kHz are, irrespective of 
the two types, marginally within the detection limits of the currently running detector of 
the first LIGO and the detection seems more feasible for the next-generation detectors such 
as LCGT and the advanced LIGO for a Galactic supernova. It is true that the GWs in 
the low frequency domains mentioned above are relatively difficult to detect due to seismic 
noises, but a recently proposed future space interf erometers like Fabry-Perot type DECIGO i s 



designed to be sensitive in the frequency regimes ( iKawamura et al.ll2006l : iKudoh et al.ll2006l ). 
The sensitivity curve of the detector is plotted with the black line in Figure [3 Our results 
suggest that these low-frequency signals, if observed, could be one important messenger of 
the increase-type waveforms that are likely to be associated with MHD explosions exceeding 

IQ51 gj,g_ 



5. Summary and Discussion 

By performing a series of two-dimensional SRMHD simulations, we studied signatures 
of GWs in the MHD-driven core-collapse supernovae. In order to extract the gravitational 
waveforms, we presented a stress formula including contributions both from magnetic fields 
and special relativistic corrections. By changing the precollapse magnetic fields and initial 
angular momentum distributions parametrically, we computed twelve models. As for the 
microphysics, a realistic equation of state was employed and the neutrino cooling was taken 
into account via a multifiavor neutrino leakage scheme. With these computations, we found 
that the total GW amplitudes show a monotonic increase after bounce for models with a 
strong precollapse magnetic field (lO^^G) also with a rapid rotation imposed. We pointed out 
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Fig. 7. — Gravitational-A yave spectrum for re presentative models wi th the expected d etection 
hmits of the first LIGO JAbbott et al.lboosi ). the advar iced LIGO JWeinsteinll2002h. Large- 
scale Cryogenic Gravitational wave Teles cope (LCGT) fjKuroda fc the LCGT Collaboration 



20061 ) . and Fabry- Perot type DECIGO (IKawamura et al 



2006 



Kudoh et al 



noted that /ichar is the characteristic gravitational wave strain defined in 



2000). It is 



Flanagan fc Hughes 



(119981 ). The supernova is assumed to be located at the distance of 10 kpc. 
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Table 5: Summary of GW amplitudes for all the models. The colors used for clarify are the 
same as Table [2l The first column represents the GW amplitudes when we terminated the 
simulation (100 ms after bounce). The second column, hpeak, is the peak GW amplitudes in 
the spectra. The third column, hio^, is the average amplitudes below 100 Hz. The supernova 
is assumed to be located at the distance of 10 kpc. 



that this trend stems both from the kinetic contribution of IvlHD outflows with large radial 
velocities and also from the magnetic contribution dominated by the toroidal magnetic flelds 
that predominantly trigger IvlHD explosions. For models with weaker initial magnetic flelds, 
the total GW amplitudes after bounce stay almost zero, because the contribution from the 
magnetic flelds cancels with the one from the hydrodynamic counterpart. These features 
can be clearly understood with a careful analysis on the explosion dynamics. It was pointed 
out that the GW signals with an increasing trend, possibly visible to the next-generation 
detectors for a Galactic supernova, would be associated with IvIHD explosions exceeding 10^^ 
erg. 

Although the presented simulations have utilized the leakage scheme to approximate the 
deleptonizati on, it would be mor e accurate (especially before bounce) to employ a formula 



developed by iLiebendorferl (120051 ). which was designed to flt ID Boltzmann results. Figure [8] 
shows snapshots at around 33 ms after bounce for model B12X20/30.1 in which the leakage 
scheme (left panel) or the prescription (right panel) is employed, respective ly. Note that 



G15" is taken in our simulation until bounce among the parameter sets in iLiebendorfer 



and is switched to the leakage scheme at the postbounce phase. As is shown, the 



^The inner-core mass at bounce for the employed parameter set is 0.6 Mq for our non-rotating 25 Mcr 



progen itor. This value is higher than that obtained in GR simulations (0.45-0.55 Mq) in iDimmelmeier et al. 
(|2008l ). This may be because the pseudo GR potential employed in this work underestimates the GR gravity, 
which could potentially lead to a large inner-core mass. 
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Fig. 8. — Snapshots at 33 ms after bounce for model B12X20/90.1 in which the leakage 
scheme (left panel) or the prescription (right panel) is employed, respectively. In each 
panel, density (logarithmic, left-half) and entropy (right-half) distributions are shown. The 
side length of each plot is 600x600 (km). 
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Fig. 9. — Left and middle panels show magnetic pressure (red line) vs. ram pressure (blue 
line) for model B12X20/30.1 around 32 ms after bounce with the Ye prescription (top panels) 
or the leakage scheme (bottom panels) along the polar axis (left panel) or the equatorial 
plane (middle panel). Matter pressure is shown by green line as a reference. The right 
panels show velocity profiles along the pole near after the stall of the bounce shock (red 
lines). Both in the two different deleptonization schemes, the MHD-driven explosions are 
indeed obtained. 
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Fig. 10. — Gravitational waveforms for model B12X1/31 near bounce with different grid 
points. Three lines starting from (at t — tb = -2 ms with tb being the epoch of bounce), 
correspond to models with different angular grid points (30 (green), 60(red), 120(blue)) while 
fixing the radial grid points to be 300. The bottom three lines are set to start from -100 
in the GW amplitudes (just for convenience), and they correspond to models with different 
radial grid points (250 (pink), 300 (orange), 600 (brown)) while fixing the lateral grid points 
to be 60. Note that the fiducial set employed in this work is 300(r)x60(6'). 
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shock revival also occurs for the model with the prescription (right panel). Figure [9] 
depicts the magnetic pressure (red line) vs. ram pressure (blue line) along the polar axis 
(left panel) or the equatorial plane (middle) near the rebirth of the stalled shock for the 
model with the prescription (top panels) or the leakage scheme (bottom panels). For the 
equator, the magnetic pressure is much less than the ram pressure (middle panel), while the 
magnetic pressure amplified by the field wrapping along the pole becomes as high as the 
ram pressure of the infalling material at the shock front, leading to the MHD-driven shock 
formation (see right panels). Regardless of the two different deleptoniaztion schemes, these 
important features associated with the MHD explosions are shown to be quite similar. 

Now we mention a c omparison between the o btained results and relevant MHD simula- 



tions. Model R4E1CF in IScheidegger et al.l ( 120101 ) whose precollapse rotational parameter is 



P = 0.5 % with a uniform rotation imposed, and whose initial poloidal magnetic field is set 
to 10^^ G, is close to our model B12X20/30.1. From their Figure 23, the jet propagates to ~ 
300 km along the rotational axis at around ~ 18 ms after bounce. In our counterpart model, 
the MHD-driven shock revives after around 30 ms after bounce, and it reaches to 300 km at 
around ~ 10 ms, which is equivalent to ~ 40 ms after bounce. Considering that our model 
(/3 = 0.1%) is a slower rotator than model R4E1CF (/3 = 0.5%), the delay of the sh o ck re- 



vival for our model seems reasonable. Model s20AlB5-D3M12 in ICerda-Duran et al.l (120081 ) 



whose precollapse angular velocity is 4 rad/s (the rotational parameter should be close to 
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P = 0.1 %) with uniform rotation and whose initial poloidal magnetic field is 10^^ G, is 
close to our model B12X20/90.1. From their Figure 10, the MHD jet propagates to ~ 240 
km at 51 ms after bounce. The employed EOS is the same as ours (t he Shen EOS), whil e 
the deleptonization scheme taken in their study was the formula (iLiebendorferl |2005| ). 
As mentioned abov e, the dynamics is rath er close to our corresponding model. Among the 
computed models in lBurrows et al.l ( l2007al ). model M15B12DP2A1H which has a precollapse 
angular velocity of tt rad/s (the rotational parameter should be close to (3 = 0.06%) with 
initial dipolar magnetic field of 10^^ G, is close to our model B12X20/30.1. The interval 
before the launch of the MHD shock for their model is 80 ms after bounce (e.g., their Table 
1) is much later than our model (28 ms after bounce). This may be due to the larger initial 
angular momentum (/3 = 0.1%) assumed in our study. 

Model A3B3G5-D3M13 in Obergaulinger et al. (2005) which has a rotational parameter 
of /3 = 0.9% with a differential rotation imposed (the radial cut off is 500 km) and the initial 
poloidal magnetic field is 10^^ G, is closer to our model B12X1/30.1. The MHD jet reaches 
to 500 km at around 7 ms after bounce, which is also the case of our counterpart model. As 
discussed above, our results are compatible to the ones obtained in the relevant foregoing 
results. 

The major limitation of this study is the assumption of axisy mmetry. Recently it was 



repor ted in three-dimensional (3D) MHD core-collapse simulations (jScheidegger et al.l 12008 



20101) that the fast growth of the spiral SASI hinders the efficient amplification of the toroidal 
fields, which could suppress the formation of jets rather easily realized in 2D simulations. As 
a sequel of this study, we plan to investigate the 3D effects in SRMHD. Regarding a resolution 
dependence of our results. Figure [10] indicates that our standard resolution is adequate to 
follow the evolution of the computed models. However, it is not suffici ent at all to capture 
the magneto-rotational instability (MRI, e.g., Balbus fc Hawlevl (1998)). At least 1 — 100 
times finer mesh points are required for resolving the MRI ( lObergaulinger et al.ll2009l ). which 
may require some adaptive- mesh- refinement treatment, a very important component that 
remains to be improved. If the MRI could be resolved, the increasing-type waveform could 
emerge also for models with weakly initial magnetic fields because a more efficient field 
amplification could be captured. Although the general relativistic effects were treated only 
by a very approximative way, we think that the general relativity should not drastically 
change our results qualitatively, because the central protoneutron stars will not collapse to a 
black hole during our simulation time as inferred from a simple argument of the compactness 
of the inner-core. As for the microphysics, the neutrino heating is not included in this study. 
However the inclusion of the neutrino heating may play a minor role in the waveforms, sinc e 



the timescale s befo re the neutrino-driy en explosions set in (e.g., iMarek &: Jankal ( 120091 ): 



Bruenn et al.l ( 120101 ): ISuwa et al.l ( 120091 ) and references therein) are much longer than the 
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MHD explosions. 

As one possible extension of this study, we think it interesting to study GWs from 



Miiller et al. 


( 


2004 


); 



total amplitudes become larger when we include the neutrino GWs. This is because the neu- 
trino GWs may make a positive contribution since the neutrino emissio ns from the oblately 



deformed protoneutron stars become stronger t oward the polar direction ( jJanka fc Moenchmeyer 



19891 : iKotake et al.ll2003al : IWalder et al.ll2005t lOtt et al.ll2008h . If this is really the case, the 



total GW amplitudes especially for the increasing type should be much larger, possibly mak- 
ing its detection more promising. Furthermore, the neutrino signals from MHD explosions, 
may have a sharp directional dependence through neutri no oscillations, reflect ing the as- 



pherical propagation of the shock to the stellar envelope (iKawagoe et al.ll2009l ). Taking a 
correlation analysis between the neutrino and GW signals could help to reveal the hidden 
nature of the central e ngines. In fact, several observational proposals have been m ade in 



this direction recently (Ivan Elewyck et al.l l2009l : lAso et al.l l2008l : iLeonor et al.l |2010| ) . The 



MHD-driven core-collapse supernovae, albeit rather minor among typical type II supernova 
explosions, seem to still contain a number of rich research subjects. 
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